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Abstract 

Gravitational coupling of the propagation four-vectors of matter wave functions is formulated in 
flat space-time. Coupling at the momentum level rather than at the “force-law” level greatly 
simplifies many calculations. This locally Lorentz-invariant approach (G4v) treats electromagnetic 
and gravitational coupling on an equal footing. Classical mechanics emerges from the incoherent 
aggregation of matter wave functions. The theory reproduces, to first order beyond Newton, the 
standard GR results for Gravity-Probe B, deflection of light by massive bodies, precession of orbits, 
gravitational red shift, and total gravitational-wave energy radiated by a circular binary system. 
Its predictions differ markedly from GR for the gravitational-wave radiation patterns from rotating 
massive systems, and for the LIGO antenna pattern. G4v predictions of total radiated energy from 
highly eccentric Kepler systems are slightly larger than those of similar GR treatments. A detailed 
treatment of the theory is in preparation. However the generation and detection of gravitational 
waves is exactly the same as the corresponding treatment for electromagnetic waves given in Collective 
Electrodynamics^ (hereinafter referred to simply as CE) and therefore separable from the material 
in preparation. It therefore seems advisable to make the gravitational-wave material available, since 
its predictions should be testable as data from Advanced LIGO becomes available over the next few 
years. The presentation is somewhat more detailed than would be “normal,” simply to make the 
approach clear and accessible to non-specialists. 


*©2015 Carver Mead 
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1 Motion Masses due to Gravitational Waves 


In G4v, the gravitational vector potential A is the direct measnre of how the local frame of reference is 
affected by the movement of distant matter. If we choose onr coordinate system for local physics to be 
dehned by the universe at large, the moving matter will have a very small influence on the coordinate 
system. In that global coordinate system, a free-floating mass in a uniform gravitational scalar potential 
will have a velocity 

V = cA_i_ (1) 

Let us consider two masses, mi located at ri and m 2 located at f ^2 relative to the binary source. The 
relative position of the two masses is / = 7^2 ~ G- Because of the vector potential of the source masses, the 
positions fi and ^*2 will be functions of time. Each mass will move around a little orbit, periodic in time. If 
the two masses are part of a gravitational-wave detector and the binary source is at astronomical distance, 
the spacing between the two masses will be vanishingly small compared to either individual distance. For 
that reason the direction and amplitude of the two motions will be the same, only executed at different 
times due to the difference in arrival time St = {r 2 — ri)/c of the potentials at the two masse^ 

ri{t) - ri(0) = f 2 {t + St) - f 2 {St) (2) 


If the wavelength of the gravitational wave is much longer than the distance between the two masses 


f2{t St) 

so the vector I from mi to m 2 is 


dr 

r 2 {t) + -^St = f 2 {t) +v St 


I = r 2 {t) - ri{t) ~ Iq- V St 
where Iq = f 2 {St) — ri(0) is not a function of time. 


(3) 


(4) 


Using R as the unit vector in the direction of ri or r 2 , the distance between the two masses thus 
becomes 

/ ~ Iq — V ■ I St 


cA ■ i St 
R-l 


h 

Iq — cA ■ i 
Iq-[ A-1 


(5) 


R 


■0 


The fractional change in length between the two masses thus becomes 

IR-l 


( 6 ) 


We can understand this formula intuitively in the following way: 

The hrst term arises because the velocity of each individual mass is in the direction of A. The second 
term expresses the fact that the position of each mass is the velocity times the time. Thus the difference 
in the positions will be the difference in vector potential times the difference in time of arrival. 


^Throughout this document, when a vector name is used without the vector arrow above it, it means the modulus 
(“length”) of the vector, and the vector name with a hat over it means a unit vector in its direction. 
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2 Gravitational Wave Radiation 


In G4v, propagating waves are characterized by the Green’s fnnction of the sonrce energy-momentnm 
4-vector, which takes the place of the charge-cnrrent 4-vector of electrodynamics. In other respects the 
calculations are exactly the same as their electrodynamic counterparts, which have been set forth in 
CE. In the far-held, which is all we will consider in what follows, the transverse vector potential is the 
only contributor to the wave’s interaction with other matter, leading to substantial simplihcation of the 
calculations, as detailed in Section 

r G Pi G 

~ A Vi^ A Ti 


2.1 Circular Binary Source 

We hrst consider the simplest possible source of gravitational waves—a binary system consisting of two 
masses Mi and M 2 orbiting their common center of mass with angular frequency G, as shown in Fig. 
Although simple, this model is an excellent approximation of many real binary systems. Mass Mi has 
velocity Vi at orbital radius ri and Mass M 2 has velocity V 2 at orbital radius r 2 , where Miri = M 2 r 2 - 



Figure 1: Top and side views of binary gravitational wave source. Masses Mi and M 2 are orbiting the y 
axis (axis of rotation). The velocities are Vi and V 2 -i and the radii of the orbits are ri and r 2 - Observation 
is done at a point located at a distance i? ri, in the y — z plane at an angle l from the y axis. The angle 
i is called the inclination of the binary with respect to the line of sight. The time difference 5ti — 6 t 2 is 
calculated from the distances c 6 ti and c 6 t 2 along the path to the observer. 

The 2 ; positions of the two masses are 

zi = ri cos Qt , , 

o, ® 

Z 2 = —r 2 cos i It 
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Thus the time differences 6 ti and 6 t 2 along the path to the observer between the Green’s functions of the 
two masses i^ 

sin t = —ri cosGt sin t , 

(9) 

c St 2 = —Z 2 sin L = r 2 cos fit sin i 

where the reference position { 6 t = 0) is taken at the center of the binary. We restrict this simple example 
to velocities v c. The z-directed velocities of Mi and M 2 are 


dzi ^ ^ 

Vzi = = — iZrismiZt 

ot 

dz2 

Vz2 = WTT = ^^2 sin fit 
ot 

The velocities vi and V 2 with respect to the observer are those given above multiplied by sin l. 


( 10 ) 


2.1.1 Vector Potential 


The vector potential A at an observation point at distance i? 3> ri is the sum of the Green’s functions of 
the individual masses (Eq. [^: 

The X components of the momenta are 


Pxi ~ MiVix = Mifl Ti cos fit = Mfl r cos fit 
Px 2 ~ M 2 V 2 X = —M 2 flr 2 Cosflt = —M fir cos fit 

The components of the momenta are 

Pzi ~ MiViz = —Mifl Ti sin fit = —M G r sin fit 
Pz 2 ~ M 2 V 2 Z = M 2 G r 2 sin fit = Mfl r sin fit 

and the positions are given by Eq. 

The last form is possible because MiVi = M 2 r 2 = Mr and pi = —p 2 . 

We will often use the quantities r = (ri + r2)/2 and M = 2 M 1 M 2 /{Mi + M 2 ), 
which reduce to the obvious mass and radius of the symmetric system. 

For historic reasons, most treatments use the reduced mass p = M 1 M 2 /{Mi + M 2 ) = M/2. 
We calculate the the x and 2 ; components of A separately: 


( 11 ) 


( 12 ) 


2.1.2 Ax component 

From Eq. 

Ax{R,t) ^ — 
(A 

^Because observation is done at a distant point 
the delay. 


Pxi {t - 6 ti)_ ^ J^x 2 {t - 5t2) 


(13) 


R R 

in the y — z plane, the difference in x makes a negligible contribution to 
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where we have neglected the contribution of the x location to the separation since S> ri, and have reset 
the zero of time to the arrival time at R. Using Eq. 11 

GMVt r / \ 

AJRR) ^—( cos — (5fi) — cosflit — 5 ^ 2 )) 

Rc^ \ / 


Rc^ 

GMQr 

Rc^ 


cos Qt cos fl 6 ti + sin fit sin fl 6 ti 


(14) 


— cos fit cos fl 6 t 2 — sin fit sin fl 6 t 2 
GMfl r / 

cos fit {cosflSti — cos flSt 2 ) + sin fit (sinfi^ti — sinfi5t2 


Rc^ 

All fl5t <C 1, so the cos terms may be neglected to hrst order and the sin terms approximated by their 
arguments: 


A^{R,t) 


GMflr 


sin fit (fi(5ti — fl 5 t 2 


Rc^ 

GMfl^r . f n r2 . 

--— sm lit-cos i it sm i -cos i It sm t 

Rc^ VC c 


GMfl^r 

RA 

GMff 

Rc^ 


sin fit cos fit sin l ^ri + 


(15) 


sin 2fit sin t 


where the second form uses 5ti and 5 t 2 from Eq. [^and r = (ri + C2)/2. 

The result can be more easily understood when the individual contributions are spelled out: 


A 

■^X ~ 


GM fir ^ 

■ Air- sm 2i it ■ sm t 


Rc^ 


(16) 


The negative sign merely indicates the phase of the sinusoidal signal due to our chosen origin of the 
rotation. The hrst term is the fraction of the frame of reference contributed by each mass M at distance 
R. The second term is the phase shift across the orbit—if both contributions were in phase they would 
cancel out at large distances and there would be no far-held radiation. The third term is the velocity of 
each mass, which it has to impart to the frame of reference. The fourth term is the time dependence with 
which the masses have momentum in the x direction and the hfth term is the dependence on inclination 
angle l of the binary system relative to the direction of observation. 

The x-directed held is given by 


^ GMfl^ 

Cx = —jT— ~ — COS 2\lt sin i 


dt 


Rc^ 


(17) 


The Transverse Theorem (Section]^ tells us that only the vector potential transverse to the direction of 
propagation contributes to interaction with matter. In the coordinate system we are using, the components 
of the transverse vector potential A± are Ax and A^. 


wr. X . GMfl^r'^ 

A^{R,t) = Az cos L ^ - - COs2iZt SlU t COS t 


RA 


A, ^ —hr, cos 2 fit sin i cost where hn = 


GMff 


(18) 


Rc* 


The t-directed held is given by 


^ dA, GMQ^r^ 

£, = —~ ——-— sin 2 flt sin 2 l 


dt 


RA 


(19) 
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2.1.3 Radiation Pattern Comparison with GR 


We are now in a position to make direct comparison of the radiation patterns predicted by four-vector 
gravity (G4v) and GR. 


90 1 90 1 



Figure 2; Comparison of quadrupole radiation patterns for G4v (left) and GR (right). Amplitude is 
plotted in the radial direction as a multiple of vs inclination l. Red-blue is the polarization for G4v, 
and the A+ polarization for GR. Magenta-cyan is the polarization for G4v, and the A^ polarization 
for GR. The phase of the red and magenta lobes is 180 degrees from that of the blue and cyan lobes. The 
zero of inclination is when the source rotation axis is pointing toward the observer. 


Both theories predict radiation at twice the rotation frequency of the source. Both theories have two 
polarizations in phase quadrature, with the same amplitudes ho'- 

GR : ^ (l -f cos^ i) cos2Gf 

GAv : A^ = ho (sin t cos l) cos 2Qt A^ 

The amplitudes of the two radiated polarizations are, however, quite different functions of the inclination 
angle l. Plots of the two polarization amplitudes as functions of inclination l are shown for both GR 
and G4v in Fig. The GR polarizations are “stretch-squeeze” tensors, while the G4v polarizations are 
ordinary vectors. GR predicts non-zero amplitude at all inclinations, with a maximum on the rotation 
axis (t = 0). G4v predicts zero amplitude on the rotation axis, with maximum amplitude normal to the 
rotation axis (t = 90°). This is a strong prediction of G4v; one of the few cases where its predictions differ 
from those of GR at hrst order. 


= ho cos L sin 2VLt 
. = ho sin L sin 2VLt 


( 20 ) 


2.1.4 Field Energy 

The gravitational vector potential is dimensionless. 

The gravitational field £ = — VV — dA/dt. has the units of inverse time. 

In the electromagnetic case, the electric held times the quantity of charge gives the force. In the G4v 
formulation of gravitation, the held £ times the quantity of matter Q = me gives the mass times the 
acceleration, which is the force. 

Thus far the parallel between the two works well. 
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In the electromagnetic case, the held energy per unit volume of a propagating wave is given by: 


= eo (^) (21) 

What is the gravitational equivalent of eo? 

The potential V of a single charge q charge is 


V = 




Aneor 

In the G4v formulation of gravitation, the potential is 

G mc^ G me 


V = 


c'i j- 


r 


GQ 

r 


The correspondence in coupling constants is therefore 

1 G 

^ or 


47ren 


eo 


47rG 


( 22 ) 


(23) 


(24) 


Therefore the held energy per unit volume Wy of a propagating wave is just that due to the transverse 
held 

The power per unit area S carried by the wave is just this energy density multiplied by the velocity c 

^3 


S=CWy = 


AttG 


G> + (i-h 


(26) 


2.1.5 Total Radiated Power 


As in the electric dipole case, the total power radiated P is obtained by integrating S over the sphere of 
radius R 


FromEq. 17 we have 


FromEq. 19 we have 


P= S 2ttR sin i Rdi 

Jo 

R^c^ r 


2G 


'0 


) sini di 


(27) 


^2 AG'^M^^Pr^ 2 r. 2 
~ -cos^ 2VLt sm^ L 




R2 c8 


sin^ 2Qt sin^ 2l 


(28) 


(29) 


Because (cos^ 2f2f) = (sin^2f2f) = 1/2, Eq. 27 for the average power becomes 


P 


Ac^ 


A sim i + sin^ 2l ) sin l dc = 


8 


( 30 ) 


which, using the average mass M, agrees with Eq. 43 when the eccentricity e = 0. 

















2.2 Eccentric Binary Source 


As we have seen for the circular binary, the Green’s function for the radiated gravitational vector potential 
A at an observation point at distance R ^ r is the sum of the Green’s functions of the individual masses 
(Eq.gi. Of course we must use the retarded Green’s functions of the two masses: 




Pi {t - 6ti) ^ P2{t- 50 ) 


R 


R 


G ^ 


G f^P G dp d - s G dp ^ ^ 

Rc^ \ dt y Rc? dt c RA dt 



(31) 


where d is the vector distance from mi to m 2 and 6t = 6ti — 6t2 is the delay of the signal from m 2 relative 
to that from mi in the direction of propagation s. We have neglected the contribution of the orbital 
location to the amplitude since R^ ri and recognized that p = Pi = —P2- For the hnal form we have 
made the approximation that the time delay 5t across the orbit is small compared with the period of any 
radiated harmonic. In PetersfcMathews |39] . following the derivation of Landau&LifshitzjlS], this is one of 
the assumptions of their quadrupole approximation. 

The gravitational field E is given by the time derivative of the vector potential: 


- dA G d f dp ^ ^ \ 

dt Rc^ dt\dt ^ ) (32) 

Es = E-s El = E‘^ -El 

where Eg is the component of the held in the direction of propagation and E^ is the component transverse 
to the direction of propagation. 


2.2.1 Total Radiated Power 

By the Transverse Theorem, the far-held vector potential is purely transverse, and therefore only the 
transverse component of the held E±_ contributes to the radiated power. From Eq. the instantaneous 
far-held power P per unit area propagating away from the source is: 

^ = ^ (33) 

darea IttG df2 AttG 

where the hnal form recognizes that the unit area darea = R^dG, where i7 is the solid angle in the 
direction s from the source (not to be confused with the rotational frequency G in the previous section). 
Eq. I^is the G4v equivalent of Peters & Mathews Eq. 3. 

The line of apsides passes through the foci of the two ellipses dehning the trajectories of mi and m 2 , 
and through the center of momentum of the system, and forms the 6 = 0 axis of the spherical coordinate 
system. We adopt the same coordinate system used for the circular binary, so the line of apsides is the z 
axis. The two ellipses he in the same plane, which forms the x — z plane of the coordinate system. The 
y axis of the coordinate system is normal to the plane of the ellipses and passes through the common 
focus of the two ellipses, which is the center of momentum. The coordinates of the two masses mi and m 2 
derived above are given in this coordinate system with the origin of the orbital phase angle "0 being the z 
axis. The vector f and the momenta of the two bodies all lie in the plane of the orbits, and therefore 
dp/dt lies in this plane as well. 
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Using the same nomenclature as Peters & Mathews, we take a as the semi-major axis of the elliptical 
orbit dehned by the distance between the two masses. In the Cartesian basis the Kepler orbital equation 
is: 

a(l-a 


r = 


and the angular velocity is: 


1 -|- e cos -0 
d%l) 


{sin-i/), 0, cos'll 




d = 2f 


We check the orbital period Pf, for consistency: 

r2iT 1 


Ph = 


'0 


= 27r^ 

dt 


G{mi + m2) 

in agreement with standard treatments. The momentum relations are thus 

^ df df dijj dp dp d'lp 

Expressing the unit vector in the direction of propagation s in spherical coordinates 

s = {sin 9 sin 0, cos 0, cos 9 sin 0} 


(34) 

(35) 

(36) 

(37) 

(38) 


we have obtained the elements of Eq. 32 Using for the time derivatives of a vector the operator 
{dip/dt)d/d^fj, and for the square the dot product of the vector with itself, Eq. 33 was evaluated using 
standard Mathematica vector operators, yielding an analytical expression for dP/dfl of considerably 
higher complexity than Peters&Mathews Eq. 14. The instantaneous total power P{ip) is then a function 
of ip, obtained by integration over the sphere: 




2tt gp 

I sin 6 dp d9 

n dil 


The output of this integration is much less complex because of the smaller number of variables: 


P{i/j) = Pcom 1 ^( 1 + e COS (?/^))^ + ysin^ (ip) 


whereas the corresponding Peters&Mathews Eq. 15 is: 


where, in both cases. 


/ ^2 

PpuiP) = Pcom i {l Pecos {ip)Y + — sin^ {ip) 
32 G^mi^m 2 ^{mi + m 2 ) (l + ^ cos {ip)) 


p = 

rom 


(39) 


(40) 


( 41 ) 


( 42 ) 


It is quite remarkable that these two approaches, starting with different fundamental theories of gravitation 
which predict drastically different radiation patterns, come up with total radiated power of exactly the 
same functional form, and differ only in the coefficient of the sin^ {ip) term. This difference may be 
detectable in systems of larger eccentricity than those yet studied. 

The average power output of the binary is obtained by integrating the instantaneous power over one 
orbital period P}y\ 




P{ip) ,, 32G‘^mi^m2^(mi-I-m 2 ) 

dip = 


b Jo 


dip/dt 


(1 - 


19e2 5e^ 

1 H-^- 

6 12 


(43) 
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which is to be compared to Peters&Mathews Eq. 16: 


+ m 2 ) / 73e^ 376^^^ 

5a5c5 (1 - e2)^/2 V 


(44) 


2.2.2 Time Derivative of Binary Period 


What is measured in observations of distant binary systems is not the energy loss, but the temporal change 
in the orbital period Pf, of the binary system. As the system loses energy due to gravitational radiation, 
the two masses sink deeper into their mutual potential well. The total energy E of the system is 


mim2G 

tj = - 

2a 


dE mim 2 G da 
dt 2a‘^ dt 


E dt a dt 


(45) 


Taylor and Weisberg|l6] quote the 1963 Peters and Mathews result, and also quote Wagoner [48]. where 


the period derivative is worked out for Kepler’s law. From Eq. ^ 

dPh 67r 


,1/2 


Ph = 271 


da 


And therefore 


(mi + m2) 

1 dPb 3 da 
2a dt 


Ph dt 


3 dE 


dt 2 (mi + m 2 ) 9t 
3a dE 


(46) 


mim2G dt 


96 G^ 17117712 (mi + m 2 ) 


5 


a"^ (1 


\7/2 


19e2 5e^ 


(47) 


where the hnal result, which uses (P) = —dE/dt from Eq. 43, is slightly different from the GR result 
(Eq. 44) in the coefficients of the higher order e terms in braces. 

For many astronomical sitnations, the semi-major axis is not easy to infer, whereas the period Pb is 
directly observable. So a more usefnl form of Eq. £7 is obtained by eliminating a using Eq. 

1 27r 1 


Pb \JG (mi -|- m2) 


Pb 


8/3 


Pb 


271 


27iY 

Pb) 277 {G{mi + m2)YP 


27r\ 


5/3 


27r 


(48) 


Pb) (G(mi + m2))4/3 


npon which Eq. £7 becomes 

dPb 1927r 


Pb = 


Pb = 


dt 

dPb 

dt 


-) 

cV 


5/3 


27r\ 


5/3 


17111712 


1 _j_ be 

6 “'“12 


(mi-|-m2)^^^ (1 — 


1927r /27r\ 


4 


5/3 


- 


mim2 


1 4- 4- 37 e" 

24 96 


Pb) {iTli + 7712)^^^ (1 — 


G4v 


GR 


(49) 


The quantity Pb is one of the standard post-Keplerian parameters |22j . and is slightly different in the 
nnmerator term from that obtained with GR (at hrst post-Newtonian order, IPN) as noted above. 
The term in the G4v expression is 4% larger than that of the GR expression, making the predicted 
energy loss rate for the double pulsar J0737-3039A/B (e = 0.08777) larger by only 10“^ than the GR 
prediction. Unfortunately this difference is too small to be resolved with data currently available [53], but 
may come in range as observations continue and new systems are discovered. For the Hulse-Taylor pulsar 
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PSR B1913+16 (e = 0.617) the G4v prediction is ~ 2% larger than the Peters&Mathews GR value. The 
measured period derivative for this system is about 1% higher than the Peters&Mathews valuejlT], but is 
“corrected” downward by a poorly-know acceleration relative to the solar system, thus making it agree 
with the GR prediction. At face value this result would place the G4v value outside the error bounds on 
the measured period derivative. Gaution is in order, however, not only because of the uncertainty of the 
correction, but because both calculations use the classical Kepler model for the binary, and, at least for 
G4v, a full model has yet to be developed. Such a self-consistent model is expected to have corrections of 
the order of this apparent disagreement, so the entire question awaits such an analysis. 


2.3 Lumped Pulsar Model 

We now consider the somewhat more speculative source of gravitational waves, as shown in Fig. a 
system consisting of two equal masses Mi = M 2 = M orbiting their common center of mass with angular 
frequency G, but located in planes of different y. 





Figure 3: Top (a) and side (b) views of lumped pulsar gravitational wave source. The y axis is the axis of 
rotation, and the radius of both orbits is r. The mass Mi is located ai y = a and M 2 is located ai y = —a 
Observation is done at a point located at a distance r, in the y — z plane at an angle i from the y 
axis. The angle l is the inclination of the binary with respect to the line of sight. 

(c) Magnihed third quadrant of side view showing construction for Eq. 51 


Although such a conhguration cannot exist as an isolated binary, Alan Weinstein suggested that a 
perturbation of this form may occur as a deformation of a neutron star, and therefore we should not 
ignore signals produced by such a conhguration in our searches for gravitational waves from pulsars. 

We can motivate the conjecture as follows: Pulsars, by their very nature, have a strong magnetic dipole 
whose axis is inclined with respect to the axis of rotation. The magnetic vector potential may be a 
non-negligible part of the momentum of any charged species in the pulsar. Therefore we must include the 
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magnetic contribution to the momentum as part of the source of the gravitational vector potential. As 
such, an appropriate abstraction of the pulsar as a source of gravitational radiation might be an axially 
symmetric momentum distribution (giving no gravitational radiation) plus a perturbation due to the 
magnetic momentum. The latter contribution can be modeled most simply as additional source masses 
located where the magnetic axis enters and exits the surface of the pulsar. 


Whether real pulsars have sufficient “lumps” at the entry and exit points of the magnetic axis is a question 
well beyond the scope of the present discussion. It should be pointed out, however, that the existence 
and magnitude of “lumps” of any kind on real pulsars is still an open question. For that reason we shall 
analyze the consequences of magnetically related lumps from the observational and energy-loss point of 
view without further apology. 


The positions of the two masses are shown in Fig. 

Zi = r cos VLt Vi = a 

Z 2 = —r cos Vtt 1/2 = —a 

The time difference 6t along the path to the observer is shown in the lower-right of Fig. 

c 6ti = —r cos Qt sin l — a cos l 
c5t2 = r cos Vtt sin t -|- a cos l 


( 50 ) 


( 51 ) 


where the reference position = 0) is taken at the center of the binary. We restrict this simple example 
to non-relativistic velocities. The z-directed velocities of Mi and M 2 are exactly as given in the binary 
case: 


Vzl 

Vz2 


dzi 

dZ2 


—fir sin fit 
fir sin fit 


( 52 ) 


The velocities Vi and V 2 with respect to the observer are those given above multiplied by sin t, also the 
same as the binary case. 


2.3.1 Vector Potential 


The vector potential A at an observation point at distance i? 3> r is the sum of the Green’s functions of 
the individual masses exactly as we do in electromagnetic propagation (CE 4.5): 


i 




G Pi 



rriiVi 


Ti 


(53) 


The X components of the momenta are 

Pxi ~ MiVix = Mflrcosflf 
Px 2 ~ M 2 V 2 X = —MVt r cos fit 
The components of the momenta are 


( 54 ) 


Pzi ~ MiViz = —Mflrsinflt 
Pz 2 ~ M 2 V 2 Z = Mfl r sin fit 


( 55 ) 


The time delays, as given by Eq. M are the only difference between the lumped case and the symmetric 
binary case insofar as the far-held radiation properties are concerned. 


We calculate the the x and 2 : components of A separately: 
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2.3.2 Aj, component 


From Eq. and Eq. 54 


A^{R,t) ^ — 


G ( pxi {t - 5t) px2 {t + 6t) 


+ 


where we have neglected the contribntion of the x location to the separation since r. 
Using Eq. 

Ax{R,t) ^ ^ (cos— 6t) — cos^l{t + 5t)^ 


Rc^ 

GMQr 


Rc^ 


cos Qt cos fl6t + sin Qt sin fl6t 
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Using the value of 6t from Eq. 51 
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Assuming fir <C c Eq. 58 becomes 
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(59) 


Comparing this result with that of Eq. we see that the lumped pulsar produces far-held radiation at 
the fundamental rotation frequency as well as at twice the rotation frequency. The amplitudes of the 
two frequency components depend on the angle a between the spin axis and the magnetic axis since 
a = Rq cos a and r = Rq sin a. 


(60) 


Ax{R, t) ~ ho sin a ( sin a sin 2flt sin c + 2 cos a sin fit cos i 

GMQ^ Rl 


where ho = 


RA 


The x-directed held £x is given by 

dAx 2GMfl3 r2 


£x = 


dt 


RA 


^ cos 2flt sin i -\— cos fit cos t j 


(61) 


2.3.3 At component 

The Transverse Theorem tells us that only the transverse vector potential contributes to interaction with 
matter. In the coordinate system we are using, the components of the transverse vector potential A_i_ are 
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The t-directed field is given by 

dA, 2GMQ^r^ 
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Figure 4: Comparison of radiation patterns of a lumped pulsar with “lumps” inclined 45° from its rotation 
axis. Two polarizations of the transverse vector potential are shown: A^ (left) and A^ (right). Vector 
potential amplitudes (as multiples of ho) are plotted in the radial direction vs inclination l. The zero of 
inclination is when the source rotation axis is pointing toward the observer. Blue is the amplitude, and 
red is the 2fl amplitude. 


2.3.4 Radiation Patterns 


The radiation patterns of Eq. ^ and Eq. ^ are plotted in Fig. |^for an arbitrarily chosen angle of 45° 
between the spin axis and rotation axis. We see that the 12 radiation can be large, and has a maximum 
directed along the axis of rotation, whereas the 211 radiation is always zero along the axis of rotation. 


2.3.5 Total Radiated Power 


As in the symmetric binary case, the total power radiated P is obtained by integrating Poynting’s Vector 


S from Eq. 26 over the sphere of radius R 


P = y S 271R sin i Rdi 

{ ((^x> + (^f>) Sint di 


2G 


(64) 


15 
















FromEq. 61 we have 
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where the cross terms average to zero due to orthogonality of Vtt and 2VLt. 


FromEq. 63 we have 
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Because (cos^ hit) = (sin^ hit) = (cos^ 2r2t) = (sin^ 211^) = 1/2, 
Eq. 1^ for the average power becomes 
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The first term, from the 211 radiation, is the same as Eq. 30, and reduces to it when the masses are in the 
same plane. The second term is from the H radiation, which has a quite different radiation pattern from 
that of the 211 radiation. We can evaluate the relative contributions of the two components by recognizing 
that = R^, where Ro is the radius of the pulsar, and a/RO = cos a: 
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The power radiated by the two frequency components is shown in Fig. It should be noted that the If 
radiation only makes a contribution larger than the 21f radiation when the angle between the magnetic 
axis and the spin axis is less than 40°. In this region, the pulsar is not a very efficient radiator. 


2.3.6 Discussion of the Results 

Several aspects of the present treatment deserve comment: 

All derivations have been done directly with the Green’s functions of the source elements, so we have 
never had to resort to any kind of multipole expansion. The 21f radiation carries the same energy as that 
predicted by Peters and Mathews for their quadrupole source, and we have every reason to refer to it in 
those terms as well. In this structure, the energy radiated by the If component is always less than that 
radiated at 2ff by a binary system formed from the same two perturbing masses. It seems that something 
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Figure 5: Power radiated from a lumped pulsar as a function of the angle of misalignment between rotation 
axis and magnetic axis. At large angles, radiation is primarily at the second harmonic of the rotation 
frequency fl. At angles less than 40°, radiation at the fundamental freqnency dominates. 


like the lumped pulsar is conceptually possible, so it is important to include the fundamental frequency in 
coherent searches for gravitational waves from known pulsars. Similar conclusions were reached recently 
by Bejger & Kr’olak|35]. However it may tnrn ont that no Hd radiation is possible from an isolated 
structure for some symmetry/conservation reason. The resolution of this question remains in the category 
of future work. 


3 Gravitational Wave Detection 

There are three distinct ways in which a gravitational wave of this type can interact with a detector: 

1. The transverse vector potential A± appears directly in the momentum of a mass in the detector, as 
discussed in the preceding sections. 

2. The scalar potential directly modnlates the speed of light in the detector. 

3. The vector potential appears directly as part of the vector velocity of light in the detector. 

In the long-wavelength limit, effects 2 and 3 are of higher order than effect 1 by a factor of order Ql/c. 
We therefore limit the present discnssion to the direct effect of the transverse vector potential on the 
positions of the two detector masses. 


3.1 Motion of Detector Masses 


From Eq. 16, and Eq. 18, the two components of the propagating transverse vector potential A± from a 
binary sonrce are 

. , . o.. . 1 , GUn^r^ 

A^ ^ riQ sm2i2t smt where Iiq = -^ 

Rc^ (69) 

At ho cos 2Qt sin l cos l 


In the gravitational-wave literatnre, ho is called the dimensionless strain. 


17 











As discussed in Section the vector potential is the direct measnre of how the local frame of reference 
is affected by the movement of distant matter. The fractional change 61 in length I between the two 
free-floating masses given by Eq. is 

y (70) 

where i? is a nnit vector in the direction from the sonrce to the detector. 

We can nnderstand this formnla intnitively in the following way: 

The hrst term arises because the velocity of each individual mass is in the direction of A. The second 
term expresses the fact that the position of each mass is the velocity times the time. Thns the difference 
in the positions will be the difference in vector potential times the difference in time of arrival. 


3.2 Detector Sensitivity 


In today’s world, credible gravitational-wave detectors use laser interferometers whose massive mirrors 
are suspended from highly compliant springs. This mounting arrangement is necessitated by the need 
to isolate the mirrors from vibrations of the earth. The axis of the interferometer is parallel to the 
earth’s surface, so, for small movements along the axis of light propagation, the mounts are virtually 
free-floating—they experience essentially no restoring force. Each arm of the interferometer contains a 
Fabry-Perot cavity, a single example of which is shown in Fig. 



LA 


/ 



U 


Fignre 6: Side view of Fabry-Perot cavity nsed in gravitational wave detector. The vector I is the length 
of the cavity along the axis of light propagation. 


We have already analyzed the effect of the transverse vector potential on I, the vector length of the cavity. 
The nnmber of fringes n of light of wavelength A in a cavity resonator is n = 2//A = 2lc/f. 


for 61/1 we obtain 

,n) 

So, to evaluate the fringe shift due to a gravitational wave, only two projection operations are needed 
for each polarization of A±: The projection of the transverse vector potential onto the centerline of the 
cavity, and the projection of the centerline of the cavity onto the direction of observation. 


Using Eq. 70 


The actual observatory instrument uses two cavities at right angles to each other, as shown in Fig. 

Because the interferometer takes the difference between the fringe shifts of the two arms, the sensitivity of 
the detector to the two polarizations is, from Eq. [m 


= (b ■ A,) (b ■ - (I2 ■ A,) (12 ■ r) 

S, = (b ■ A,) (b ■ R^ - {12 ■ A,) (/; ■ r) 


18 










Figure 7: Cartoon of LIGO interferometric gravitational wave detector, which has two Fabry-Perot cavities, 
labeled “light storage arm” in this drawing. The ontpnt of the photodetector is proportional to the 
difference between the phases of the two arms, and hence is sensitive to a signal that “stretches” the 
distance between the mirrors in one arm and “sqneezes” the distance between the mirrors in the other. 
From the LIGO web site [33]. 


These sensitivities are the response of the detector to a vector potential of unit amplitude in the x and l 
directions of Fig. The actual fringe shift due to a signal from a distant source is obtained by using the 
actnal vector potentials and given by Eq. 69 These two vector components of the total potential 
A are called the polarizations of the wave. The ontpnt V of the LIGO detector is the snm of the two 
polarizations, each multiplied by its respective sensitivity: 


V = Aa;Sx + ~ HqISx sin 2Vtt sin i + cos 2Vtt sin l cos l 


(73) 


3.3 G4v VS GR Antenna Patterns 

The sensitivity of a LIGO detector for the two polarizations of G4v and GR are compared in Fig. As 
the Earth rotates, these antenna patterns sweep over the sky and determine the amplitude of the received 
signals, including sign changes. 

It can be seen by inspection that, if the position of the sonrce in the sky is known, the response of a 
detector to the G4v signal will be very different from that of a GR signal. This topic is analyzed in detail 
by Isi et.al. m- 

3.4 General Nature of the Interaction 

It is interesting that the G4v patterns are almost complementary to the GR patterns. In particular, the 
G4v patterns have no sensitivity on the 2 ; axis (directly above the detector), where both GR patterns 
have maximum sensitivity. We can understand the basic nature of the G4v gravitational wave interaction 
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Figure 8: Antenna patterns for the two polarizations of G4v (left) and GR (right). The observatory 
antenna arms lie along the x and y axes. The photodetector output due to a gravitational wave of unit 
is plotted as the radial distance from the origin to the surface as a function of the direction of incoming 
wave propagation. 


with matter in the following manner: The change in length 51 between the two masses separated by a 


vector distance I due to a gravitational vector potential A is, from Eq. |70 

51^ - (A-1]{R-1 


(74) 


where i? is a unit vector in the direction of wave propagation. 

The hrst term arises because the velocity of each individual mass is in the direction of A. In the far-held, 
the only components of the vector potential that survive are perpendicular to the direction of propagation. 
Thns, in terms of its effect on matter, the wave is a “shear wave”, where all elements of matter in a plane 
perpendicular to the direction of propagation move in elliptical orbits within that plane. Thus elements of 
matter will remain at constant spacing in any plane perpendicular to the direction of propagation. 


The second term expresses the fact that the position of each mass is the velocity times the time, so the 
time for a particnlar element of matter to be displaced is different due to the difference in propagation time 
along the direction of propagation. Because the motion is perpendicular to the direction of propagation, 
elements of matter will remain at constant spacing (to hrst order) along the direction of propagation. 

Thus the difference in the positions will be maximum at 45° to the direction of propagation. So, at any 
given time and distance along the path of propagation, elements of matter will be closer together along 
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one 45° direction, and farther apart in the orthogonal 45° direction. Thus, within any plane at 45° to the 
direction of propagation, matter will be alternatively “stretched” along one axis and ” squeezed” along 
the orthogonal axis. This “quadrupole” pattern is similar to that predicted by GR, but rotated from 
the axis of propagation by 45°. Unless the absolute position and polarization phase of the gravitational 
source is known, the two theories cannot be distinguished by their instantaneous effect on “freely floating” 
matter. 


Because observatories like LIGO have two arms perpendicular to each other, whose displacements are 
subtracted, the antenna patterns shown in Section 3^ are exactly the same as the “stretch-squeeze” 
of matter along the corresponding orthogonal axes. Hence these antenna patterns have four lobes, as 
appropriate for a quadrupole gravitational wave. 


As we have shown, a propagating transverse vector potential produces “stretch-squeeze” motion in free- 
floating matter much as a tensor wave would, but in planes at 45° to the direction of propagation. The 
foregoing development is completely parallel to the corresponding derivation for electromagnetic waves 
given in [CE]. It has been known for many years that electromagnetic waves, which are vector in nature, 
propagate perfectly well from quadrupole sources [43]. 

Binary pulsar systems are, in principle, capable of exquisitely precise tests of any theory of gravitation. 
Such tests rely on equally precise models of the binary system using the theory being tested. Such models 
are presently being evolved using GR, and we hope that future versions will include G4v as well. 


3.4.1 Field Quantization 

It should be noted that there is a tendency to identify quadrupole interaction with tensor propagation. 
That identity is only appropriate for theories with quantized radiation helds. In G4v and GE, the 
electromagnetic and gravitational helds are four-vector in character, and neither is quantized. In fact, 
there is not a single shred of evidence for quantization of either the electromagnetic or the gravitational 
radiation held. The Gasimir ehect is often pointed to as “proving” the existence of virtual pairs, and 
therefore of held quantization. It was shown by Jaffe|42] that the “Gasimir forces can be computed without 
reference to zero point energies. They are relativistic, quantum forces between charges and currents.” 
Virtual pairs thus have the role of simplifying certain calculations rather than providing evidence for held 
quantization. It is easy to see that virtual pairs, and therefore held quantization of the pro forma kind 
that is routinely taught in contemporary physics classes, cannot be physically real. If they were, they 
would contribute a factor of 10^^° more to the energy density of the universe than can be estimated from 
any method of observation. See[S2] for a thoughtful review of this well-known cosmological constant 
problem. Thus the hctitious energy of vacuum huctuations is seen as a deep conceptual problem in 
prevailing theory rather than a numerical or “hne tuning” problem, as it is often portrayed. 

In CE and G4v, the degrees of freedom of matter are not separate from the degrees of freedom of the held. 
Quantization occurs in the region where held waves interact with matter waves; not in the held in the 
absence of matter. In this view of the relationship of matter and held, there are no vacuum huctuations, 
and the whole cosmological constant problem does not arise. 
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4 Sky Searches 


As the Earth rotates, the antenna patterns of Fig. sweep over the sky and modnlate the amplitude of 
the received signals, inclnding sign changes. Coherent searches for pulsar signals heterodyne the received 
LIGO signal against a sine and cosine reference locked to the pulsar radio signal (see GarchingjS^ and 
LIG0[3H] papers for details). The assumed amplitude response patterns are then used as “matched hlter” 
templates on the heterodyned signal to bring it up out of the noise. Often these searches integrate the 
signal, multiplied by the template, for a year or more. 

To interpret a signal from a particular source at a particular observatory, we need the sensitivity as a 
function of the position of the source and the phase of Earth rotation. 


4.1 Detector Coordinates 

The detector has a latitude $ (from equator), east longitude A, and i/'d is the angle between the y arm of 
the interferometer and the local northerly direction. For example Technical Note LIGO-T980044-08 gives 
the directions of the x and y arms for the LIGO Hanford and Livingston observatories. The orientation for 
Hanford corresponds to a right-hand rotation from y = north, x = east by 'ipd = 180 — 54 = 126°, and that 
for Livingston to a right-hand rotation from y = north, x = east by "0^ = 180-1- 17.7165 = 197.7165°. 

4.2 Source Coordinates 

The source coordinates are its declination a (from equator), its right ascension /?, its polarization angle 'ip 
and its inclination, shown as l in Fig. The time dependence has been combined in the local Sky Angle 
7] of the source from the observatory. The sky angle advances from to east to west with time. A source 
with larger right ascension (3 appears farther to the east, and thus has a smaller sky angle when viewed 
from a given point on Earth. The same source appears farther to the west and therefore has a larger sky 
angle when viewed from a point of larger east longitude A. Therefore ri{t) = \ — (3 + Get, where Gg is the 
earth rotation rate and t is the local sidereal time. 

A rotating source such as the binary shown in Fig. has two additional angles: 

1. The Inclination Angle l is the angle from the line of sight to the rotation axis of the source. A 
source with l = 90° has its rotation axis perpendicular to the line of sight. A source with t = 0° has 
its rotation axis pointed at the Earth. 

2. The Polarization Angle pj is the left-hand rotation angle about the line of sight from Earth to 
the source by which the source rotation axis appears to be rotated from the plane of celestial North. 

4.3 Coordinate Transformations 

To transform the detector coordinates into the source coordinate system, we hrst rotate around the z 
(vertical) axis by an angle —ipd so the x axis is pointing due east, and the y axis due north. We then 
rotate around the x (east) axis by the observatory latitude <F to align the north-pointing y axis with the 
Earth’s polar axis. Next we rotate around the y (now polar) axis by an angle —y (negative because a 
right-hand rotation goes from west to east) to take out the Earth’s rotation and place the 2 ; axis in the 
meridian of the star. Then we rotate around the x (east) axis by an angle —a to point the z axis at the 
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star. Then we rotate around the 2 ; axis by an angle —ip to align the y axis with polarization angle of the 
star. Finally we rotate about the y axis by vr to transform into the source frame, where the z axis is 
identihed as the propagation vector R, pointing at the Earth, and the y axis is in the plane of the source 
rotation axis, as shown in Fig. |T| After these rotations, the result is a vector containing the projection of 
the arm unit vector onto A^, and R directions in the source frame of reference. 


h • Ax 


■ 1 ■ 

h ■ A, 

= Ry{'K)RR-'ilj)Rx{-a)Ry{-rii)RxR)Rz{,-'ipd) 
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1 

.h-R. 


0 


(75) 


(76) 


where the brackets indicate vectors, the i?’s are rotation matrices about the respective axes by the angle 
indicated, and evaluation proceeds from right to left. 


Note that these R rotation matrices rotate the coordinate system (alias rotation), not the vector in a hxed 
coordinate system (alibi rotation), and are thus the transpose of the usual alibi rotation matrices. 


Because the interferometer takes the difference between the fringe shifts of the two arms, the sensitivity of 
the detector to the two polarizations is, from Eq. 


= (/l ■ A,) ■ R^ - (/2 ■ A,) (/2 ■ r) 

s, = {li ■ A^I (/; ■ idj - (I 2 ■ A,) (/2 ■ r) 


(77) 


These sensitivities are the response of the detector to a vector potential of unit amplitude in the x and l 
directions of Fig. The actual fringe shift due to a signal from a distant source is obtained by using the 
actual vector potentials and A^, given by Eq. ^ These two vector components of the total potential 
A are called the polarizations of the wave. The output V of the LIGO detector is the sum of the two 
polarizations, each multiplied by its respective sensitivity: 


V = AxSx + AiS't ~ ho (^Sx sin 2Qt sin t cos 2Qt sin l cos t j 


(78) 


In general the actual data will be a series of discrete samples Vi recorded at times R. Since the source 
inclination l is not a function of time over the periods involved in the search, we write Eq. 73 for V) as a 
function of R 

(79) 


Vi = ho sin 2^1R + bi cos 2VLT^ 


where a* = SxiRi) sin l and bi = S^iRi) sin l cos l. 

If Uj = bi, Eq. 79 would represent a circularly polarized wave. Because, as shown in Fig. bi < ai 
always, the G4v signal is elliptically polarized. 


5 Continuous-Wave Searches 

Any spinning mass distribution that is not uniformly distributed around its rotation axis will radiate 
gravitational waves of the kind described, and will therefore be potentially detectable by a LIGO-type 
detector. The most promising candidates are those that are the closest to us, have the most “lopsidedness”. 
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and about which we have some knowledge of the rotation rate A class of objects that meet these 
criteria are the pulsars that have resulted from recent cataclysmic astronomical events. The reception of 
electromagnetic pulses from these objects makes it possible to construct a time base that is presumably 
related to the spin frequency, and therefore to the gravitational-wave frequency kl. 

The approach to detecting gravitational waves from these objects was pioneered by the Garching group |37] 
and has been rehned and extended by members of the LIGO and Virgo collaboration. See in particular 
Matt Pitkin’s Thesis [3B] and references therein. 

The technique starts with a heterodyne (multiplication) of the original data with both a sine and a cosine 
function, with frequency cu chosen to he at twice the pulsar rotation rate. These two heterodyned series 
are then lowpass hltered (for reasons of antialiasing) before resampling. We will refer to the real functions 
resulting from the heterodyning of the original data V as the cosine and sine quadratures: 


Cj = Vi cos uiTi = hnitti sin 2GTj -|- hi cos 2GTj 


= ho ^ ( sin (2G -t- Ci;)Tj -t- sin (2G — u:)T^ 


+ hQ—[ cos (2G -|- oj)Ti -|- cos (2G — ijj)Ti 


Si = Vi sin wTj = ho sin 2GTj -|- hi cos 2GTj 

= ho-^ ( cos (2G — u:)Ti — cos (2G -|- a;)Tj 


-|- ho — sin (2k2 -|- ijj)Ti — sin (2k2 — uj)Ti 
The low-pass hlter step removes the sum-frequency terms, leaving 

Ci = ho Y sin (2G — uj)Ti + cos (2G — uj)Ti 

Si = ho-^ cos (2k2 — u)Ti — sin (2k2 — uj)Ti 


(80) 


(81) 


The radio astronomy community has developed exquisitely precise methods for synchronizing time 
coordinates. Applying those techniques to the pulsar problem has made possible the creation of precise 
sine and cosine waveforms that are locked to the electromagnetic “beeps” received from the pulsar. 


We thus have at our disposal during the heterodyne step of Eq. a time coordinate such that wT, = 
2GTj — 0, where the phase 0 between the electromagnetic “beeps” and the gravitational wave is not 
known. 


5.1 Search Techniques 

The Ah sample output of the heterodyne process contains a cosine part c* and a sine part Sj 


Ci • J, I 


Si tti , hi . 


( 82 ) 


where a* and hi, given by Eq. 79, are periodic functions of the time Ti with period 1 sidereal day. 


The most straightforward technique for recovering a signal from random noise is the matched filter, in 
which we multiply each signal sample by a template representing the sensitivity for a unit-amplitude 


24 




version of the expected signal at the time Tj of that sample, sum the results over all data samples available, 
and divide by the number of samples. The result is an estimate of the amplitude of the expected signal 
actually present in the data. The two signals we have are c* and Sj, which we can combine as an analytic 
function ^ 

° tti sin (j) + bi cos 0 + f (uj cos (j) — hi sin ( 


Vi= [d + lSi] = — 


If we set 
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Eq. 1^ become^ 
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(83) 
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and our matched-filter template Fj is 


F,- = 


H _ \/Q? + H 




( 86 ) 


An estimate of the amplitude of the signal can be formed by taking the dot product of Vi and T), treated 
as vectors. A source with the opposite rotation will reverse the sign of the imaginary part of the signal. 
For that reason, a match was computed for both right-hand (RH) and left-hand (LH) rotation: 
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h(<f) = (3}(!.,)3}(F,) - a(t>.)9(F,) 
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LH 


(87) 


Because the phase 0 is not known, the calculation of Eq. ^ is carried out for a number of values for 0, 
typically 36. Experimentally, a plot of h as a function of 0 is an extremely smooth sine wave, even for 
data that is pure noise. For that reason we take the maximum of such a plot minus the minimum as our 
estimate h^. of the amplitude hg: 


he = -^max(h(0)) — min(h(0)) 


( 88 ) 


5.2 Search Data Preparation 

The process began with heterodyned data prepared with exquisite care by Matt Pitkin. The largest 
outliers were trimmed as described in Section At this point the standard deviation of the complex data 
for HI was 1.59 x 10“^^. Before attempting a search for any actual signal, it was important to establish 
the background statistical properties of the data as it is reflected in accidental correlations between the 
noise and the templates. For this purpose we require a version of the received signal in which there is no 
Crab signal. 

^This result could have been obtained immediately by treating both the radiated signal and the antenna patterns as 
analytic functions from the beginning. 
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This Catch-22 situation is solved in the following way: The complex data are re-heterodyned by multiplying 
each complex data sample by where t is the sample time and h is a frequency well within the passband 
of the heterodyned data. This operation does not change the amplitude, but approximately orthogonalizes 
the signal with respect to the template. A data set of this type was prepared using frequencies from .008 
Hz below the Crab frequency to .008 Hz above, spaced by 1.16 x 10“® Hz, giving a total of 1.4 x 10® 
frequencies. Each frequency was re-heterodyned to base band and matched to both right-hand and 
left-hand rotation G4v templates. The results are shown in Fig. 



Figure 9: Cumulative distribution of detected amplitudes for 1,400,000 searches of HI S5 data at frequencies 
spaced evenly above and below the Crab “beep” frequency. The form of this plot was chosen to emphasize 
the tail of the distribution. 


5.3 Signal Insertions 


A set of 100 G4v Crab signals was prepared with amplitudes randomly distributed from zero to 3% of the 
standard divination of the data, and phases randomly distributed from 0 to 27r. A subset of the original 
1.4 X 10® data sets were augmented with 100 additional sets, to each of which was added a different one 
of the synthesized signals. The results are shown in the left plot in Fig. By keeping track of the 


amplitude of each insertion, we can compare it with the amplitude actually recovered by the search of the 
data set containing that insertion, as shown in the right plot in Fig. 


These two plots give us two methods to estimate the probability that a detected amplitude could have 
occurred by chance: For example, the often-used rule-of-thumb that any point greater than 5.326724 a 
has only one chance in ten million of arising by chance would place the detection level at 1.635 x 10“^® 
using the right plot. The same criterion would place the level of detection at 1.615 x 10“^® using the 
expression for the red line in the left plot. These detection thresholds are just above 1% of the standard 
deviation of the original input signal. 
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Figure 10; Left; Cumulative distribution of detected amplitudes for 2400 searches. 2300 of these data 
sets are a subset of those shown in Fig. |^and contained no added signal. 100 additional sets contained 
synthesized signals randomly distribnted in amplitnde and phase. Right; Comparison of the amplitnde 
recovered in the search vs the amplitnde actnally inserted for the 100 additional data sets 

A much more comprehensive search process, including searches sensitive to both GR and G4v signals 
and method for distingnishing between them is given in Isi ef.a/.[M]- There it is shown that “For some 
combinations of detectors, sonrces, and signal strengths, G4v signals are invisible to GR templates and 
vice versa.” Thus the imminent detection of gravitational waves by Advanced LIGO and other highly 
sensitive observatories will give us a powerful method for comparing theories of gravity in this new dynamic 
domain. 


6 Power Line Contamination 

The LIGO Crab paper|l9] contains the following assertion; “Althongh 2u [59.56 Hz] is close to the 60 Hz 
power line frequency, it is sufficiently far away that the searches are relatively unaffected by nonstationary 
components of the power line noise.” It is the pnrpose of this section to examine the issne of power-line 
contamination in more detail. The issue is complex for a number of reasons. Although the long-term 
average freqnency of the U.S. power grid is held at 60 Hz to high accnracy, instantaneons deviations can 
occnr in local areas dne to sndden increases or decreases in load. Power-line noise therefore has a strong 
daily variation due to the life and work patterns of our human society. These time dependencies are 
evident in the load cnrves kept by electric power providers, as illnstrated by Fig. [TTj 

Factories are started up every morning. Lights come on every evening. Induction motors run at a rotation 
rate lower than the synchronous speed (60 or 30 Hz) by an amount dependent on their load. Under no-load 
conditions a high-quality induction motor can run arbitrarily close to synchronous rate, progressively 
slowing as load is increased. Florescent lights generate large amonnts of noise near 60 and 120 Hz, with 
considerable bandwidth due to the inconsistent hring point of the plasma as the applied AG voltage 
waveform changes sign. All of these phenomena make predicting the noise properties near line freqnency 
problematic. 
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Figure 11; Load as a function of time for a region in the southwestern Unted States. The left curve is 
for one week-day in Jannary and the right cnrve is for one week-day in Angnst . In both plots the daily 
maximum is more than twice the daily minimum. Notice the remarkable difference between the total 
power usage and the daily pattern in the two plots. This huge daily variation and seasonality of power 
nsage is typical of many areas in the United States. 

6.1 Statistical Properties of the S5 Crab Data 

The hfth LIGO science run (S5), started on 2005 November 4 and ended on 2007 October 1. During 
this period the detectors (the 4 km and 2 km detectors at LIGO Hanford Observatory, HI and H2, and 
the 4 km detector at the LIGO Livingston Observatory, LI) were at their design sensitivities and had 
duty factors of 78% for HI, 79% for H2, and 66% for LI. No gravitational-wave signal was detected 
from the Grab Pulsar at a level above 1% of the standard deviation of the norm of the heterodyned 
data[l9]. Thus the outliers must be from causes other than the Crab, and any correlation with the Crab 
gravitational-wave signal must be purely accidental. For that reason we only consider the two most 
sensitive detectors HI and LI. 




Figure 12: Cumulatve distribution functions of HI (left) and LI (right) heterodyned data norms from the 
S5 Crab search. In both plots the red and dark blue points are the fnll data set and the light bine points 
are for the same data with the largest outliers (red points) removed. The reduced sets contain 90% of the 
heterodyned data, corresponding to a dnty factor of 70.2% for HI and 59.4% for LI. 




































The cumulative distribution functions of the heterodyned signal norm from these two detectors are shown 
in the Fig. This plot format was chosen to best portray the behavior of the “long tail” of outliers. It 


is clear that the distribution is unusual, with a lot of weight carried by the outliers, which could not have 
been due to the Crab signal. Alan Weinstein suggested excluding the worst offender^ and re-analyzing 
the data to see if, perhaps, the outliers, normally associated with transient events, were also affecting the 
long-term spectral properties of the search. 


6.2 Spectral Properties of the S5 Crab Data 

Any gravitational-wave signal from the Crab Pulsar will be manifest as a periodic variation of the received 
signal with period 1 per sidereal day, with a strong component at 2 per sidereal day. Unfortunately, that 
period is within 1 part in 365 of 1 per solar day. 

If power-line contamination were strictly periodic with solar-day period, we could separate it cleanly from 
the sidereal-day period gravitational-wave signal by using one year’s worth of data. However we are not 
that fortunate for several reasons. The S5 run was not quite a year in duration (201 days of data for HI, 
and 158 days of data for LI). Power-line load patterns vary a great deal on a yearly timescale. Summer 
entails a large air-conditioning load that is not present in the winter. Daylight savings time shifts the 
phase of some human activity and not that of others, etc., etc. These yearly variations, evident in Fig. [T^ 
modulate contamination signals, and thus create sidebands shifted from the solar day period by exactly 
the same spacing as the sidereal day signal 

Because the Earth rotates around its axis in the same sense that it rotates around the Sun, the number of 
solar days in a year is one less than the number of sidereal days in a year. Because of the trigonometric 
identity 2 cos A cos H = cos {A + B) -(- cos {A — B), any signal periodic at one per solar day and ampli¬ 
tude modulated at one per year will have sidebands at ±l/(days per year) from the one per solar day 
fundamental. The upper of these sidebands will therefore coincide with the sidereal day frequency. 

The spectral estimation process for data of this kind is problematic because the data are not uniformly 
sampled. The nominal sampling rate after heterodyning and low-pass hltering is 1 per minute. However 
there are many gaps in the data where no samples exist. Standard Fourier transform techniques cannot 
be applied directly under these conditions. We were fortunate to contact Dick Lyon of Google just 
after he had updated the Wikipedia page on Least Squares Spectral Analysis. He suggested using the 
Lomb-Scargle algorithm, which is specihcally optimized for nonuniformly sampled data. Matlab code and 
source reference are given in Section 

We begin by analyzing the Lomb-Scargle spectrum of each observatory in the neighborhood of the first 
two harmonics of 1 per day, as shown in the following plots. The most obvious difference between the two 
observatories is that all the LI harmonics are much larger than those from HI. The second thing to note 
is that all the spectra are changed in fundamental ways when 10% of the data, made up of the largest 
outliers, have been removed. 


“Trimming” is the term often used in the statistical literature for this procedure. 
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Figure 13: Spectrum of the norm of HI (Left) and LI (Right) heterodyned data in the neighborhood of 
one per day. The red curves correspond to the full datasets and the blue curve is from the 90% subsets 
shown in light blue in Fig. 12 The main peak is reduced considerably when 10% of the data, made up 


of the largest outliers, are removed. The peak around one per sidereal day, critical for any signihcant 
detection, is reduced by a factor of 2 in the HI data and by nearly a factor of 10 in the LI data by this 
minor trimming process. 




Figure 14: Spectrum of the norm of HI (Left) and LI (Right) heterodyned data in the neighborhood of 
two per day. The red curves correspond to the full datasets and the blue curves are from the 90% subsets 
shown in light blue in Fig. 12 The main peak is rednced considerably when 10% of the data, made up 


of the largest outliers, are removed. The peak around one per sidereal day, critical for any signihcant 
detection, is reduced by a factor of 2 in the HI data and by nearly a factor of 10 in the LI data by this 
minor trimming process. The sidebands out to the sidereal second harmonic are greatly rednced in the 
data sets that exclude the largest outliers. Note that the LI spectrum in this region is about an order of 
magnitude larger than that of HI. 
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6.3 Discussion 


There are many phenomena that have a daily period modulated on a yearly basis, the power line load 
being only one of them. However the narrow (1/60 Hz) bandwidth of the Crab heterodyned data near 60 
Hz makes power-related noise particularly suspect. If this supposition is taken as a working hypothesis, it 
might be possible to hnd correlation between some of the LIGO environmental monitors and the outliers 
identihed in the present study. On the other hand, it is entirely possible that the Crab heterodyne is itself 
the most sensitive indicator of disturbances in its narrow band of frequencies. There are an uncountable 
number of ways the power line can couple into any experiment. The LIGO interferometer, being the most 
sensitive instrument ever devised by the human race, is vastly more prone to contamination than any 
other. The consequence of this situation is that no environmental monitor is likely to be anywhere near as 
sensitive to this contamination as the interferometer itself. 


6.4 Hanford vs Livingston 

The integrity of the power grid near Hanford is particularly high. The facility is near the center of one of 
the most extensive hydroelectric power developments in the worldjSn]- Not only do the nearby power 
plants have a large excess generating capacity, but, more important in the present context, they provide a 
massive “spinning reserve”—the term for synchronous, electrically coupled angular momentum that lends 
stability to the phase of the power grid voltage in the face of load transients. Perhaps nowhere else in the 
world is the power grid environment so favorable for a near-60 Hz search. These circumstances lead us to 
interpret the remarkable superiority of the Hanford Crab data as being, most probably, the result of its 
superior power-line integrity. 

No gravitational-wave signal was detected from the Crab Pulsar at a level above 1% of the standard 
deviation of the norm of the heterodyned data[l9]. Thus the outliers must be from sources other than 
the Crab, and any correlation with the Crab gravitational-wave signal must be purely accidental. This 
being the case, the transient outliers in the raw heterodyned data can be discarded as described here with 
negligible impact on the statistical signihcance of the result due to “snooping the data” isn. Calculating 
the exact revision to the p-values determined in the search due to this data cut is left to the experts. 

We conclude that discarding at least 10% of the data points, corresponding to the largest outliers in the 
heterodyned Crab data, is a logical and reasonable step in the Crab Pulsar search process, and that any 
effect on the integrity of the result due to this step is likely to be be negligible. 
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7 Appendix - Transverse Theorem 


We present this theorem for electromagnetic radiation. 

In G4v, the relations for weak gravitational radiation are identical with two changes: 

1. The charge-cnrrent-density 4-vector is replaced by the energy-momentnm 4-vector. 

2. The conpling constant is —Gjc? rather than l/47reo. 

The Green’s Fnnction integral forms for the retarded scalar and vector potentials are: 


V{t) = 


1 f p{t- r/c) 


dvren 


dvol 




J{t — r/c) 


dvol 


( 89 ) 


where r is the distance from the sonrce point to the point at which V and A are measured. 

We consider a compact charge distribution that is changing with time, and inquire about the general 
properties of the far-held radiation from that charge distribution. By far-held we mean that the distance 
r from any point in the source to the point of observation is larger than any dimension of the source by a 
large enough factor that, for any two points in the source, 1/r is constant to within the accuracy required. 
We choose our coordinate system so that the vector from the source to the point of observation is in the x 
direction. We divide the source into thin slices with planes of constant x, such that each slice of thickness 
6x contains charge q{x,t)- We choose the origin x = 0 such that its distance to the point of measurement 


is R. From Eq. 89 the potential V contributed by the charge g(x, t) in a single slice located at x at time t 
will be delayed by the propagation time {R — x)/c 




47reoi? 


( 90 ) 


Therefore the potential Vi contributed by the charge qi (x, t) in a single slice located at x is 

1 


Vi(i') = 


dvreni? 


9i(«) 


(91) 


where t = t' — R/c R x/c. The potential at the same observation point from the immediately adjacent 
slice (2) located aX x + 5x is 

'^ 2 (i') = D <? 2 (t + Sx/c) (92) 


dvreni? 


So the total potential V at the point R will be 


V = Vi + V2 = 


dvreoi? 


+ q 2 (t + 5x/c] 


dV 


/ N / N dqo 

A'KCqR \ c ot 

1 cdqi dq 2 Sx d‘^q 2 

AneoR \ dt dt c dR 


(93) 


If we consider these two adjacent slices in isolation, the charge on the hrst will be decreased and that on 
the second will be increased by a current I crossing the boundary between them in the -|-x direction. 


By conservation of charge 


dq2 _ _ j 


dt 


dt 


(94) 
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which, with Eq. 93 gives 


dV 


6x d‘^q 2 


AirenRc dR 


( 95 ) 


The potential V is propagating in the +x direction at velocity c. Its fnnctional form is therefore 
V = f(t — R/c). The derivatives are thus related by 


M 


1 W 

c df 


Sx d‘^q2 jj^oSx d‘^q2 
attr or 


AneoRR dR 


( 96 ) 


where we have used the fact that = 1/c^- 


By Eq. 89 and Eq. 94, the vector potential at point R will be 


A/D+'\ ^0 fJ^oSx dq2 

dA nodx d‘^q2 
dt AttR dR 


(97) 


From Eq. 
vanishes: 


96 


and From Eq. we conclude that the total electric held S in the direction of propagation 

(98) 


dA _ dV 

We can make up any compact charge distribution out of a superposition of pairs of slices like we have just 


£ = - VV - ^ = 0 
dt 


analyzed. Because the Green’s Functions (Eq. 89) are linear, if the longitudinal electric held is zero for 


each pair of slices individually, it will be zero for the sum as well. 

We have thus arrived at a very important result for an electromagnetic wave propagating in free space far 
from its source: 


Transverse Theorem: The electromagnetic wave propagating in free space due to a compact source has 
the following property: 

The time derivative of the longitudinal vector potential just cancels the spatial gradient of the scalar 
potential; therfore the electric field in the direction of propagation Zero. 


The far-held radiation from any source approaches a plane wave arbitrarily closely at large distances from 
its source. The scalar potential V in a plane wave is constant within any plane normal to the direction of 
propagation. It follows from the transverse theorem that the only attribute of the wave that can affect 
charges is the transverse vector potential A±, and the only non-zero component of the electric held is 
T_l = —dA±/dt. Plane electromagnetic waves are thus said to be Transverse Waves. From the Green’s 
Function (Eq. 89) it follows that we can hnd the components of the vector potential that interact with 
matter by considering only the components of source current normal to the direction of propagation. This 
deep result enables an enormous simplihcation in calculating any far-held radiation pattern. The special 
case of this result for an electric dipole source is shown in GE 4.16. 


This conclusion is with regard to propagating wave solutions. It obviously does not apply to the static 
electric held of a charge distribution with a net charge that is not varying with time. 
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8 Appendix - Lomb-Scargle Algorithm 


Matlab code for the Lomb-Scargle algorithm as given in; 

Ice Ages and Astronomical Causes: 

Data, spectral analysis and mechanisms 

By Richard A. Muller, Gordon J. MacDonald 

\protect\vrule widthOpt\protect\href{http://books.google.com/books?id=P8ideTkMQisC}{http://bO' 
+Iomb+scargle&hl=en#v=onepage&q=spectraiyo201ombyo20scargle&f=f alse 

function power=lomb(t,y,freq) 
y Lomb-Scargle periodogram 
nfreq=length(freq); 
fmax=freq(nfreq); 
fmin=freq(l); 

power=zeros(nfreq,1); 

f4pi=freq*4.*pi; 

pi2=2*pi; 

n=length(y); 
cosarg=zeros(n,1); 
sinarg=zeros(n,1); 
argu=zeros(n,1); 
var=cov(y); Zvarlance 
y subtract mean 
yn=y-mean(y); 
for fi=l:nfreq; 

sinsum=sum(sin(f4pi(fi)*t)); 

cossum=sum(cos(f4pi(fi)*t)); 

tau=atan2(sinsum,cossum); 

argu=pi2*freq(fi)*(t-tau); 

cosarg=cos(argu); 

cosnorm=sum(cosarg.*cosarg); 

cfi=sum(yn.*cosarg); 

sinarg=sin(argu); 

sfi=sum(yn.*sinarg); 

sinnorm=sum(sinarg.*sinarg); 

power(fi)=(cfi*cfi/cosnorm+sfi*sfi/sinnorm); 

ypor normalized output 

Zpower(fi)=(cfi*cfi/cosnorm+sfi*sfi/sinnorm)/(2*var); 

end 

The plots in this document were made without normalizing the output 
in order to compare the noise magnitude at HI and LI observatories. 
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